set terminal epslatex color standalone lw 1 header '\renewcommand{\normalsize}{\scriptsize}  \usepackage{txfonts} \newcommand{\dd}{\mathrm{d}}'
set output"M_vir_tmc3.tex"
power(x)=a*(x**p)    
fit [0.05:40] power(x) 'tmc.txt' using 9:10 via 'power.par'
yequal(x)=x
set ylabel offset 4,0
set key box
set key width -23
set key height 0.4
set xrange [0.07:40]
set yrange [1:100]
set logscale xy
unset key
set xlabel '$M_{LTE}\ (M_\odot)$'
set ylabel '$M_{vir}\ (M_\odot)$'
plot"tmc.txt" using 9:10 w p pt 19 title "Data" , power(x) with line lt 1  title \
'$\frac{M_{vir}}{M_{\odot}}=(22.9\pm 2.7)\left(\frac{M_{LTE}}{M_{\odot}}\right)^{0.39\pm 0.04},\ R^2=0.77$',\
yequal(x) with line lt 4 title '$M_{LTE}=M_{vir}$'

set output"M_j_tmc.tex"
line(x)= k*x+b
power(x)=a*(x**p)    
fit [0:150] power(x) 'tmc.txt' using 9:11 via 'power.par'
set xlabel '$M_{LTE}\ (M_\odot)$'
set ylabel '$M_{J}\ (M_\odot)$'
plot"tmc.txt" using 9:11 w p pt 19 title "Data" , power(x) with line lt 1  title \
 '$\frac{M_{J}}{M_{\odot}}=(17.0\pm 2.8)\left(\frac{M_{LTE}}{M_{\odot}}\right)^{0.35\pm 0.06},\ R^2=0.75$',\
yequal(x) with line lt 4 title '$M_{LTE}=M_{J}$'
set terminal wxt enhanced
set output

set terminal epslatex color standalone lw 1 header '\renewcommand{\normalsize}{\scriptsize}  \usepackage{txfonts} \newcommand{\dd}{\mathrm{d}}'
set output"M_vir_cmc.tex"
power(x)=a*(x**p)    
fit [0.05:150] power(x) 'cmc.txt' using 9:10 via 'power.par'
yequal(x)=x
set ylabel offset 4,0
set key box
set key width -23
set key height 0.4
set xrange [1:150]
set yrange [10:400]
set logscale xy
unset key
set xlabel '$M_{LTE}\ (M_\odot)$'
set ylabel '$M_{vir}\ (M_\odot)$'
plot"cmc.txt" using 9:10 w p pt 19 title "Data" , power(x) with line lt 1  title \
'$\frac{M_{vir}}{M_{\odot}}=(22.9\pm 2.7)\left(\frac{M_{LTE}}{M_{\odot}}\right)^{0.39\pm 0.04},\ R^2=0.77$',\
yequal(x) with line lt 4 title '$M_{LTE}=M_{vir}$'

set output"M_j_cmc.tex"
line(x)= k*x+b
power(x)=a*(x**p)  
set xrange [1:150]
set yrange [10:400]  
fit [1:150] power(x) 'cmc.txt' using 9:11 via 'power.par'
set xlabel '$M_{LTE}\ (M_\odot)$'
set ylabel '$M_{J}\ (M_\odot)$'
plot"cmc.txt" using 9:11 w p pt 19 title "Data" , power(x) with line lt 1  title \
 '$\frac{M_{J}}{M_{\odot}}=(17.0\pm 2.8)\left(\frac{M_{LTE}}{M_{\odot}}\right)^{0.35\pm 0.06},\ R^2=0.75$',\
yequal(x) with line lt 4 title '$M_{LTE}=M_{J}$'
set terminal wxt enhanced
set output


